Overview of MHz air shower radio experiments and results 



Benoit Revenu 

Ecole des Mines de Nantes - Universite de Nantes - CNRS/IN2P3 
4 rue Alfred Kastler, BP20722, 44307 Nantes CEDEX03, FRANCE 

Abstract. In this paper, I present a review of the main results obtained in the last 10 years in the field of radio-detection of 
cosmic-ray air showers in the MHz range. All results from all experiments cannot be reported here so that I will focus on 
the results more than on the experiments themselves. Modem experiments started in 2003 with CODALEMA and LOPES. In 
2006, small-size autonomous prototypes setup were installed at the Pierre Auger Observatory site, to help the design of the 
Auger Engineering Radio Array (AERA). We will discuss the principal aspects of the radio data analysis and the determination 
of the primary cosmic ray characteristics: the arrival direction, the lateral distribution of the electric field, the correlation with 
the primary energy, the emission mechanisms and the sensitivity to the composition of the cosmic rays. 
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1. INTRODUCTION 

In 1971, Allan [1] gave a review of the latest results, at 
that time, concerning the detection of extensive air show- 
ers initiated by high-energy cosmic rays. The initial mo- 
tivations for studying the radio signal was the possibil- 
ity to build an array of widely spaced receivers to de- 
tect extensive air showers at the highest energies (above 
10 EeV), in a cheap way. It was also proposed to use the 
radio signal in conjonction with a particle array to get 
additional information on the longitudinal profile, pro- 
viding valuable constraints on the nature of the primary 
cosmic ray. Using the data available at that time, AUan 
concluded that the main signal should be polarized in the 
direction of the vector v x B, where v is the direction 
of the shower axis and B the direction of the geomag- 
netic field. The mechanism responsible for this specific 
polarization is the geomagnetic contribution due to the 
Lorentz force acting on each secondary charged parti- 
cles, in particular the electrons and positrons. It was pre- 
dicted that the charge excess contribution should be less 
important but still detectable, mainly for incoming direc- 
tions close to the direction of the geomagnetic field. It is 
also stated that the electric field amplitude extrapolated 
on the shower axis, over the range 32-55 MHz, is pro- 
portional to the primary energy. Finally, for a given event 
observed by several detectors, the electric field amplitude 
decreases exponentially with the axis distance. The equa- 
tion proposed to fully describe the electric field Ey ob- 
served by a single detector at a distance R of the shower 
axis is: 



in jiV m^ MHz^ , where Ep is the primary energy, 
a the angle between the shower axis and the geomag- 
netic field, 6 the zenith angle, R the distance between 
the detector and the shower axis and Rq the attenuation 
length of the electric field. The lateral distribution func- 
tion (LDF) proposed by Allan depends on the axis dis- 
tance leading to an azimuthal invariance of the electric 
field with respect to the shower axis. This formula was 
used as a starting point at the beginning of the years 2000 
for the design of the CODALEMA [2] and LOPES [3] 
experiments. These first modem experiments were trig- 
gered by an array of particle detectors in order to search 
a posteriori for a radio counterpart to a shower candi- 
date. The next challenge consisted in setting up fully au- 
tonomous and independent radio detectors to test the pos- 
sibility to build a full stand-alone radio array. A first type 
of such radio array was installed in 2006 at the center of 
the surface detector (SD) of the Pierre Auger Observa- 
tory, in the scope of the Auger Engineering Radio Array 
(AERA), experiment which finishes its phase 1 in 2012. 
AERA benefits from the SD and FD (fluorescence detec- 
tor) reconstructions, providing high-quality super-hybrid 
events. The LOFAR [4] experiment also detects the radio 
emission of air showers. The complementary particle de- 
tector, LORA [5], helps in triggering and identifying the 
cosmic rays detected by LOFAR. 

We present in this overview the updated results of the 
field. We will discuss the lateral distribution function, the 
correlation with the primary energy and we will insist 
on the current situation concerning the emission mecha- 
nisms that can be determined through the polarization of 
the total electric field. 



2. DETECTION OF THE ELECTRIC 
FIELD EMITTED BY EAS 

The secondary electrons and positrons created during the 
development of the air shower form a pancake-shape par- 
ticle front moving at the speed of light. The thickness of 
the shower front is of the order of 1 m on the shower axis, 
up to ^ 10 m far from the axis. These particles suffer 
a systematic opposite drift caused by the Lorentz force 
due to the Earth's magnetic field generating a coherent 
emission of electromagnetic waves in the 1-500 MHz 
range. The amplitude of the resulting macroscopic cur- 
rent and its variation depends on the number of charges 
particles of the shower and therefore, can provide infor- 
mation on the longitudinal profile. The electric field pro- 
duced this way is expected to have a polarization fol- 
lowing the direction of the Lorentz force, v x B. This 
mechanism is known as the geomagnetic contribution as 
discussed in the introduction. Another mechanism is due 
to the variation of the excess of electrons due to the an- 
nihilation of positrons and knock-on electrons, known 
as the charge-excess contribution. This corresponds to 
the Askaryan [6] effect in the air, leading to a radially- 
polarized electric field, also in the MHz domain. It is 
therefore possible to disentangle these two contributions 
through the different polarizations of the associated elec- 
tric fields. Various simulation codes are available, a re- 
view is given in this conference, see [7]. This electric 
field is usually detected by systems composed of anten- 
nas, precise time tagging system (GPS for instance) and 
the signal is digitized by fast ADCs in order to be able 
to compute the Fourier spectrum up to some hundreds of 
MHz. The electric field is filtered in a frequency band 
above the AM (~ 20 MHz) and below the FM (80 MHz). 
Following the antenna used, the filtering can be done in 
more restricted bands accordingly to the frequency re- 
sponse of the device used. In general, the electric field is 
measured in, at least, two horizontal orthogonal polariza- 
tions; this permits to reconstruct the polarization angle 
in the horizontal plane and to check for the underlying 
emission mechanism. 

The two first modem radio experiments, CO- 
DALEMA and LOPES, led to many progresses in 
the understanding of the emission of the electric field 
by EAS. These experiments detect the EAS initiated by 
cosmic rays with an energy between 10^^ and 10^* eV 
and follow the same principle: they use a particle array 
which triggers the radio array. The particle array for 
CODALEMA is an array of 17 scintillators and in the 
case of LOPES, the particle array is the KASCADE- 
Grande experiment. LOPES is installed close to the 
city of Karlsruhe which is an electromagnetically noisy 
environment. The solution adopted by LOPES is to use 
interferometric technics, using as an input the high- 



quahty information provided by KASCADE-Grande. A 
digital beam-forming is applied to inter-phase the time 
series of all LOPES antennas using the arrival direction 
given by KASCADE-Grande. Then, a cross-calibration 
procedure is used to estimate the global amplitude of 
the electric field. For high signal-to-noise ratio events, 
the electric field can be estimated for each antenna, this 
permits to compute the LDF of these events. On the 
contrary, the CODALEMA experiment is installed in 
the radio observatory of Nangay which is a protected 
area in the sense that the electromagnetic emissions of 
the neighbourhood is controlled. This permits to use the 
data of the radio array of CODALEMA independently 
of the data of the particle array, once the EAS is clearly 
identified and validated by the particle array. Each radio 
detector provides the electric field as a function of time 
with a high sampling rate (typically 500 Ms/s or 1 Gs/s) 
so that it's possible to measure the time of the maximum 
of the signal. 

3. ARRIVAL DIRECTION 

From the measurements of the electric field by several 
detectors, it is possible, in the very same way than with a 
particle array, to reconstruct the incoming direction of 
the EAS. Provided the ground coordinates of at least 
three non-aligned detectors and the time of transit of the 
electromagnetic wave of the shower in these detectors 
has been measured, we can triangulate to estimate and 
<j), the zenith and azimuth angles of the shower, respec- 
tively. For higher multiplicity events, we can compute 
the radius of curvature of the electromagnetic wave as- 
suming a spherical front. The time resolution of GPS re- 
ceivers allow to reach an angular resolution of the order 
of a fraction of degree, even for low-multiplicity events. 
Artificial sources are commonly used to time-calibrate 
the radio detectors. For instance, the detection of air- 
plane transients represents a gold mine because it's pos- 
sible to inter-calibrate the detectors, to study the antenna 
lobe sensitivity and also to compute the angular resolu- 
tion. It has been proven that this angular resolution can 
reach 0.5°, even for a small number of detectors spaced 
by 140 m only [8]. Octocopter flights have been used for 
instance in LOFAR, LOPES and AERA. The arrival di- 
rection estimation is important at two levels: it permits to 
identify - and suppress - anthropic events (mainly com- 
ing from the horizon, apart for airplanes) and to identify 
cosmic ray events by a comparison with the arrival di- 
rection given by the particle array. If the directions agree 
within a certain angle (typically less than 20°) and the 
time of the shower is the same within a certain time win- 
dow (less than 100 ns for instance), then the event is seen 
by both arrays. 



4, LATERAL DISTRIBUTION 
FUNCTION 

The LDF describes the electric field amplitude as a func- 
tion of the distance to the shower axis. The starting point 
for the function describing the LDF is an exponential, 
as first proposed by Allan. Most of the radio-detected 
events are well described by an exponential profile (see 
Eq. 1). Nevertheless, a non-negligible fraction of events 
(of the order of 20%) present a flattening for detectors 
close to the shower axis, in particular for inclined show- 
ers, as reported by the LOPES collaboration in [9]. The 
same observation holds for the CODALEMA data and 
the LOFAR collaborations. A very nice example of such 
event is given in [4] . This event was detected by 5 LO- 
FAR stations corresponding to more than 200 indepen- 
dent measurements of the electric field at different axis 
distances; the radio pulse power as a function of axis dis- 
tance is presented in Figure 1 . The flattening of the pro- 
file is very clear for antennas located at less than 100 m 
from the shower axis. 



property is predicted by all the simulation codes and has 
also been observed in the CODALEMA data [13]. More 
complicated 2D LDFs should be used in the future anal- 
yses but at the time of this conference, no analysis using 
2D LDF has been presented. The minimization of the 
based on the model s{d) = SQexp{—d /do) gives the core 
position through the axis distance d =\nx CA| where n 
is the shower axis and CA is the vector between the core 
position and the position of antenna A. The minimization 
also provides the attenuation length do and the on-axis 
signal So- The on-axis signal can be defined as the actual 
electric field if the amplitudes s{d) have properly been 
deconvoluted for the antenna response (arrival direction 
and frequency gain). This deconvolution is mandatory in 
order to be able to study the correlation between electric 
field and energy of the primary cosmic ray. Using simula- 
tions, it has been demonstrated that the shower-to-shower 
fluctuations are minimized when considering the electric 
field at a distance of ~ 1 10 m of the shower axis. It could 
be therefore interesting to consider this value instead of 
the on-axis electric field [14]. 
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FIGURE 1. Air shower detected by the LORA/LOFAR ex- 
periments. The flattening of the profile is interpreted as an ef- 
fect of the air refractive index leading to a Cerenkov ring close 
(less than --^ 100 m) to the shower axis. 



The flattening of the electric field profiles can be under- 
stood as the effect of the air refractive index. At the sea 
level, n ^ 1 .0003 and this value decreases with increas- 
ing altitude; the consideration of realistic values of the 
air refractive index in the simulations permits to repro- 
duce events with a flat profile close to the shower axis 
which is the consequence of a Cerenkov ring (see [7], 
[10], [11], [12]). One should be cautious when using the 
usual ID exponential profile, in particular when consid- 
ering the electric field measurements close to the shower 
axis. Moreover, the electric field should not exhibit az- 
imuthal invariance with respect to the shower axis. This 



5. CORRELATION WITH THE 
PRIMARY ENERGY 

Hybrid arrays usually use the particle detector array for 
the energy estimation. As stated before, the LOPES ex- 
periment relies on the KASCADE-Grande reconstruc- 
tion, the CODALEMA experiment uses the scintillator 
array and the AERA experiment depends on the Auger 
SD and/or FD reconstruction. Allan proposed a linear re- 
lation between the on-axis electric field and the primary 
energy. This relation is confirmed by the LOPES 30 data, 
using the pulse height in the east-west (EW) polarization: 
gEw °^ £0'^5±o.04 (ggg |-j5]) ]y[Qj.g recently, the AERA 

experiment reported the correlation between the electric 
field and the primary energy determined by the Auger 
SD and FD. The signal is recorded in the three polariza- 
tions EW, north-south (NS) and vertical (V). The AERA 
radio detectors have been fully calibrated and the mea- 
surements are deconvoluted to estimate the 3D electric 
field vector [16]. Then the Hilbert envelope is computed 
for the three directions and the total signal strength is de- 
fined as the maximum of the 3D Hilbert envelope [17]. 
The electric field is finally interpolated at 110 m of the 
shower axis, where the energy resolution is maximum, 
and the correlation with the primary energy is compati- 
ble with a linear relation. The AERA energy correlation 
plot will be published in a forthcoming paper. The Pierre 
Auger collaboration also reported a positive correlation 
at 99.99% CL between the on-axis electric field and the 
primary energy, obtained with the (pre- AERA) RAuger 
setup, as can be seen in Figure 2 extracted from [18]. 
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FIGURE 2. Correlation of the extrapolated on-axis electric 
field with the primary energy estimated by the Auger SD. The 
Pearson correlation coefficient is O.SI^q^^ at 95% CL, taking 
into account all sort of systematic and statistical errors. 
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It has been reported during this conference the quasi- 
linear correlation of the CODALEMA data with the pri- 
mary energy using raw data (not deconvoluted for the an- 
tenna response) [19]. Due to the different estimators used 
in the community (either the on-axis electric field or its 
value at ^ 110 m), one should be cautious when com- 
paring the slopes of the correlations reported by several 
experiments. 



6. EMISSION MECHANISMS AND 
POLARIZATION 

The geomagnetic contribution to the total electric field 
is dominant, as stated in the 1960s and confirmed with 
much more statistics and better data some years ago 
by CODALEMA [20], [21] and LOPES [15]. In the 
southern hemisphere, the RAuger prototype reported 
that the arrival directions of the detected air showers 
were in good agreement with a v x B effect. Using 
the AERA data, we studied the correlation in the hor- 
izontal plane of the expected polarization angle 0g of 
the geomagnetically-induced electric field with the mea- 
sured polarization angle (j)p. The correlation is excellent, 
as can be seen in Figure 3 (extracted from [22]). 



FIGURE 3. Polarization angle in the horizontal plane of the 
detected electric field as a function of the expected polariza- 
tion angle in case of a pure geomagnetic contribution. Green 
triangles correspond to regular atmospheric conditions. Yellow 
circles correspond to events for which the atmospheric mon- 
itoring was not working. Red triangles corresponds to events 
detected during thunderstorms. 
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FIGURE 4. Polarization pattern at the ground level for a ver- 
tical shower in a geomagnetic field oriented along the NS axis. 
The left figure corresponds to the geomagnetic contribution and 
the right figure, to the charge-excess contribution. 



In the last year, there have been two different approaches 
to detect such an additional contribution, compatible 
with the charge-excess contribution. 



As discussed in the introduction, the electrons in ex- 
cess relatively to the positrons implies a net electric field 
whose polarization pattern is radial in the shower trans- 
verse plane. At first order, the total electric field is ori- 
ented following V X B but it is possible to search for sec- 
ond order effect, that can be revealed either by a change 
in the polarization angle or by an enhanced electric field 
amplitude due to constructive interferences. Figure 4 
presents the polarization patterns for a vertical shower 
in a geomagnetic field with a NS component. 



6.1. Evidence for a radial contribution 
through polarization angle studies 

The polarization angle of the total electric field is es- 
timated using the measurements along the EW direction 
(x-axis) and the NS direction (y-axis). To detect an elec- 
tric field contribution incompatible with the geomagnetic 
mechanism (i.e. having an orientation not aligned along 
V X B), we first construct an observable R characteriz- 
ing the deviation from a pure geomagnetic electric field. 



for each radio detector For that, we first rotate the co- 
ordinate system so that the new jc-axis, y! , is aligned 
along the direction of v x B projected in the horizon- 
tal plane and the new y-axis, /, is perpendicular to x'. 
R quantifies the relative signal strength in the direction 
perpendicular to the geomagnetic expectation. By con- 
struction, a pure geomagnetic electric field leads to 7? = 
and an electric field with no geomagnetic component has 
R~±\. For a given set of detected showers, we can 
compute /?data for the data (we used the data of a pre- 
AERA prototype) and compute the same factor 7?simu for 
the same simulated showers (same geometry and detec- 
tor array). The simulations used here are MGMR [23] 
and REAS3 [24] and were run using realistic showers 
(having an excess of electrons) and for showers forced to 
have the same number of electrons and positrons. The re- 
sults showed that the correlation is much stronger when 
including the charge-excess mechanism in the shower. 
Figure 5 presents the correlation between /?data and /?simu 
with the charge-excess taken into account for the simula- 
tion code MGMR for example. 
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FIGURE 5. Correlation between -Rdata and i?simu for ths case 
of the code MGMR taking into account the charge-excess 
mechanism in the computation of the simulated electric field. 
The correlation is clear and much better than for the case where 
the charge-excess is not included in the simulation code. 

Quantitatively, the reduced decreases from 6.4 (4.7) 
to 2.7 (3.0) when taking into account the charge-excess 
in the simulated showers with REAS3 (MGMR). The 
conclusion of these results is that the electric field emit- 
ted by the showers detected by this pre-AERA proto- 
type is in better agreement with the simulations when 
the charge-excess mechanism is included in the simu- 
lation codes. The relative influence of both prototypes 
depends strongly on the shower arrival direction and on 
the relative position of the observer with respect to the 
shower core. More details on this analysis can be found 
in [25]. The same analysis is currently performed using 
the AERA data and more refined simulation codes. 



6.2. Evidence for a radial contribution 
through core shift 

The CODALEMA collaboration reported in [13] the 
measurement of a shift toward the east of the shower core 
position estimated using the radio data (the radio core) 
with respect to the shower core estimated using the par- 
ticle array data (the particle core). The particle core po- 
sition is obtained using the Nishimura-Kamata-Greisen 
(NKG) lateral distribution (see [20] for the details). The 
radio core position is obtained by fitting an exponen- 
tial function of the type sid) = 5oexp(— J/iio) where the 
core coordinates are hidden in the axis distance d. The 
shift to the east is characterized by the quantity = 
■*^c-.radio " -^cpaiticie • This shift is interpreted as the result of 
the superposition of the geomagnetic and charge-excess 
electric field components, as presented in Figure 4. For 
this example of a vertical shower, an observer located at 
the east of the shower core will observe a constructive 
superposition of the two components of the electric field, 
contrary to an observer located at the west of the shower 
core. The radio core will then be reconstructed, in this 
example, to the east of the actual shower core, defined as 
the intersection of the shower axis with the ground. The 
shift strongly depends on the incoming direction through 
the geomagnetic component. It is therefore interesting to 
study the shift as a function of the EW component of 
V X B (because CODALEMA measures only the EW po- 
larization). To check the influence of the charge-excess, 
simulations were run (using the code SELFAS2 [26]) 
with and without the charge-excess mechanism. With no 
charge-excess, no core shift is observed. Therefore, the 
core shift is an evidence for the charge-excess contri- 
bution. For the selected CODALEMA data set, the de- 
pendence of A( with (v X B)ew is presented in Figure 6 
and compared with expectations from simulations with 
SELFAS2, ran on showers having the same characteris- 
tics than the selected showers detected by CODALEMA. 



In the last two years, there have been two evidences of 
a non-geomagnetic contribution to the total electric field. 
The additional contribution appears to be radially polar- 
ized in the transverse shower plane, as it is the case for 
the charge-excess contribution. More data on a wider en- 
ergy range will be needed to confirm if this contribution 
can be associated unambiguously to the charge-excess 
contribution (Askaryn affect in the air). The AERA ex- 
periment should be able to study in great details this con- 
tribution. 
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FIGURE 6. Shift between the radio core and the particle 
core on the EW axis as a function of (v x B)ew The data are 
represented by the circles (in blue for internals events with a 
very reliable reconstruction, in red for external events with a 
less reliable particle array reconstruction). The shaded zone is 
the ± 1 CT region determined by the simulation with the charge- 
excess contribution included. No shift is observed on simulated 
data with no charge-excess. 



7. SENSITIVITY OF THE RADIO 
SIGNAL TO THE NATURE OF THE 
PRIMARY 

The accurate determination shower by shower of the 
composition of the cosmic rays at ultra-high energies is 
the next challenge to be taken up. Since some years, us- 
ing simulations, the radio signal is expected to be corre- 
lated to the nature of the primary cosmic ray (see [27], 
[14]). From simulations, it is possible to estimate the 
value of the atmospheric depth of maximum develop- 
ment of showers X^ax. using the LDF of radio profiles. 
The LOPES collaboration report an uncertainty on X^ax 
of the order of 150 g cm^^. This large value is domi- 
nated by the very noisy environment in Karlsruhe and 
one could expect (using simulations) a much better reso- 
lution — around 30 g cm^^ — in quiet sites. For compar- 
ison, the Auger FD has a resolution of 20 g cm^^ [28]. 
We can expect, by combining data from radio and the 
SD, to reach a resolution on close to that of the 
FD. Another possible method to constrain the nature of 
the primary has been presented during this conference 
in [29]. The principle of this method is to study the slope 
(spectral index) of the frequency spectra between 40 and 
60 MHz. Simulations show that this spectral index de- 
pends not only on the shower geometry but also on the 
nature of the primary cosmic ray, at a level of 10% ac- 
cording to MGMR. The dependence on the shower ge- 
ometry has been confirmed on the AERA data and the 
possibily to constrain the composition with this method 
is still under study. 

Experimentally, the LOPES collaboration re- 
ported [30] a first evidence, at a level of 3.7(7, of the 



sensitivity of the radio signal to the nature of the primary 
cosmic ray through the longitudinal profile development. 
The slope of the radio lateral distribution (using an 
exponential profile with an estimation of the amplitude 
at an axis distance of 100 m) is correlated to the mean 
muon pseudorapidity which is in turn correlated to the 
shower development: the pseudorapidity is small (large) 
when the muons are produced at low (high) atmospheric 
heights. The slope is defined as the factor a in the lateral 
distribution function: s{d) = ^loo exp(— fl(ii— 100 m)), 
where s is the electric field amplitude expressed in 
m-i MHz-' in the band [40-80] MHz and d is 
the axis distance. Figure 7 shows this correlation using a 
selection of 59 events. 
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FIGURE 7. Mean muon pseudorapidity as a function of the 
lateral slope of the 59 selected LOPES radio events. 



8. CONCLUSION 

The most recent and important result of the field is the 
evidence for a secondary electric field, radially polarized 
with respect to the shower transverse plane. This electric 
field could well be the contribution of electrons in ex- 
cess in the shower, known as the Askaryan effect in the 
air. The confirmation of such observation and its quan- 
tification would be a major advance in the understanding 
of the emission processes of electric field in air showers. 
The radio signal, supposed to be correlated to the longi- 
tudinal profile of the shower, is very promising in the es- 
timation of the composition of ultra-high energy cosmic 
rays. The AERA experiment, located at the north-west of 
the Auger SD, it installed at an ideal site because showers 
can be detected by many detectors: the Infill Auger SD, 
the regular Auger SD, the regular Auger FD, AMIGA 
(dedicated to the muonic contents of the shower, [25]) 
and HEAT [25] which permits to detect the fluorescence 
light at higher elevation angles corresponding to lower 



energy showers. All these detectors will study in great 
details the showers in the energy range 
providing high-quality data to study the transition from 
a galactic to an extra-galactic origin of cosmic rays. The 
phase 2 of AERA, providing a total of 161 stations spread 
over 20 km'^ is starting in March 2013. 
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